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ANALYSIS OF TEMPERATURE-DEPENDENT NEUTRON TRANSMISSION 
AND SELF-INDICATION MEASUREMENTS ON TANTALUM 
AT 2-keV NEUTRON ENERGY 
by Thor T. Semler 
Lewis Research Center 
SUMMARY 

The method of pseudo-resonance cross sections is used to analyze an experiment 
which measured the temperature-dependent change in neutron transmission and self- 
indication ratio in the unresolved keV neutron energy region. The self-indication ex- 
periments consist of energy- and temperature-dependent measurements of transmission 
samples of tantalum followed directly by a capture sample of tantalum in a gamma-ray 
detector. The self-indication ratio is then the ratio of the number of captures in the 
capture sample when the transmission sample is in the beam to the number of captures 
when it is out of the beam . 

The observed transmission and self- indication ratio have been converted into effec- 
tive total cross sections for the various samples. These show variations with energy 
due to local variations in the strength of resonances. At 2 keV the energy resolution of 
the tantalum experiment is 192. 0 eV and contains on the average about 45 resonance 
levels. A general technique based on pseudo-resonance parameters is described that 
is used to extract local strength functions and a consistent set of resonance parameters. 
For the 2-keV experiment a 0. 7-meV value of the reduced scattering width "T® and a 
potential cross section Op Qt of 9. 5 barns fitted all six of the temperature -dependent 
transmission experiments and four out of six temperature -dependent self- indication 
experiments. 

Finally, a set of resonance parameters which adequately represent the temperature- 
dependent behavior of the measured effective cross sections at 2 keV is derived. This 
set is suitable for reactor calculations in this energy region. Similar analyses for the 
measurements in other unresolved energy regions are necessary to obtain local reso- 
nance parameters for improved reactor Doppler calculations. 

An evaluation of another published experiment using a scandium-filtered beam of 
neutrons upon tantalum is presented. The precision of this experiment was unfortunate- 
ly degraded by the presence of neutrons streaming down the duct containing a scandium 
filter. However, by using the measured scandium-filtered beam spectrum and recently 
measured cross sections for tantalum, some evidence is adduced for the adequacy of 
the pseudo-resonance approach for the analysis of this experiment as well. 



INTRODUCTION 


The Doppler effect is a measure of the change of effective neutron cross section 
with temperature and is of major importance to the design and safety considerations for 
fast reactors. A strong negative Doppler effect implies an increase in neutron capture 
with increasing reactor temperature and is desirable to limit prompt supercritical ex- 
cursions. Thus, an accurate estimate of the Doppler effect of reactor materials is 
needed for fast- reactor design. 

Temperature-dependent continuous neutron cross sections, in the energy region 
where individual resonances are unresolved in cross section measurements, may be 
used to calculate neutron Doppler effects in such regions. If it is assumed that the 
statistical distributions of resonance parameters obtained from analysis of the resolved 
resonance levels measured at neutron energies, from 0 to about 1 keV, are representa- 
tive of the statistical distributions at the higher unresolved energies, it is possible to 
sample from these distributions and generate "pseudo" resonance parameters (refs. 1 
to 4). These pseudo-resonance parameters may then be used to compute energy- and 
temperature -dependent cross sections from which have been evaluated the Doppler effect 
(ref. 1), the average capture cross sections (ref. 2), or the effects of overlapping fis- 
sion resonances (ref. 3). Rarely have such calculations been checked against experi- 
ments which measure the variation of the cross section as a function of temperature and 
energy in the unresolved resonance region. 

Two experiments have been performed which measured the change in neutron trans- 
mission of tantalum as a function of sample temperature (refs. 5 to 7). It was shown by 
Byoun, Block, and Semler (ref. 5, hereinafter referred to as RPI 1971) that the results 
of neutron transmission and self- indication experiments on tantalum at room tempera- 
ture could be analyzed by the method of pseudo-resonance chains. However, an earlier 
attempt to fit the experimentally measured change in effective total cross section at 
2 keV due to temperature by pseudo-resonance parameters was unsuccessful (refs. 6 
and 7, hereinafter referred to as INC 1968). 

This report demonstrates a use of the method of pseudo- resonance chains to elicit 
a set of resonance parameters which reproduce the temperature-dependent RPI 1971 
measurements at 2 keV and are suitable for reactor calculations in this region. As the 
average effective total cross section of tantalum shows considerable structure in the 
unresolved region (see fig. 5), it is clear that there are significant differences in the 
local values of the strength function. Thus, a general method which allows extraction 
of the local strength functions and consistent resonance parameters is used. The region 
at 2 keV has been chosen to illustrate this technique since there appears to be a signifi- 
cant "window" in the total cross section which should provide a severe test of the meth- 
od and since there exist other temperature-dependent data in this energy region (ref. 6). 
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This method, presently applied to tantalum, should be generally applicable to energy 
regions where the p-wave and inelastic cross sections are small or are accurately 
known. A discussion of the average resonance parameters measured for tantalum and 
used as constraints in the analysis is included. A brief critique of the INC 1968 experi- 
ment is also included. 


METHOD OF PSEUDO-RESONANCE CHAINS 

The assumption has been generally made that the statistical distributions found in 
the resolved resonance region are valid at higher neutron energies (refs. 1, 2, and 4). 
These distributions of resonance parameters are used to generate a set of resonance 
parameters in the unresolved resonance region. The code DAISY (ref. 8) allows one to 
generate chains of pseudo-resonance parameters in the unresolved resonance region 
which preserve the correct average values as well as accurately represent the correct 
distribution functions. 


Wigner Level Spacing Distribution 

The level spacing between two neutron resonances, of like J and l, follows a 
probability distribution of the form 


p(x) = - x exp 
2 



x — 0 


> 


p(x) = 0 


x < OJ 


( 1 ) 


where J is the spin quantum number of the compound nucleus formed by the target 
nucleus and the neutron, l is the angular momentum of the bombarding neutron with 
respect to the target nucleus, and x is a normalized dimensionless level spacing or 


x 



( 2 ) 


where S is the nearest neighbor level spacing in eV and Dj is the J- state average 
level spacing in eV. 

When the Wigner level spacing distribution is written as equation (1), the expecta- 
tion of x is 1. In order to sample X, a level spacing, from this normalized Wigner 
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distribution, we form the functional mapping (eq. (3)) 


R 


n 



( 3 ) 


where R n is a random number uniformly distributed on the interval (0, 1) which gives 


R 


n 



dx 


(4) 


Solving for X and substituting R fi for (1 - R n ), we obtain 

X = y- In (R n ) (5) 

We may now easily gain the correct J-state level spacing by multiplying X by Dj, the 
J- state average level spacing. 


Level Width Distribution of Porter and Thomas 

0 2 

The variation of the reduced scattering widths T is found to follow a x proba- 
bility distribution function with one degree of freedom for s-wave {1=0) resonances, 
while two degrees of freedom have been suggested for p-wave {l = 1) resonances; that is, 


r 


0 

n 


= e¥ 2 


u } T n 
L n 


( 6 ) 


where Eq is the resonance energy in eV, v ^ is the penetrability factor, and T n is 
the scattering width associated with the Breit-Wigner formula. The distribution of the 
dimensionless reduced scattering widths 


y 



(7) 


is found to be of the form 
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(8) 


p m (y) 


m 



my\(m/2)“l 
— ^ exp 

i 2 ) 



where the number of degrees of freedom m is equal to 1 for s-wave resonances and 

2 for p-wave resonances. This analysis will be restricted to s-wave resonances since 

the strength function and penetrability for p-wave neutrons at 2 keV for tantalum are 

quite small. This is the distribution of Porter and Thomas for scattering widths. 

o 

It can be shown that the problem of sampling an element from a x distribution 
with m degrees of freedom can be reduced to the problem of sampling m random 
values drawn from a normal distribution of mean 0 and variance 1. That is, if the se- 
quence 


Y v Y 2 ,Y 3 ,Y 4 ,. . .,Y r ,. . . 

is composed of values selected from a normal distribution of unit variance and mean 0, 
the sequence of transformed random variables 


y2 tt 2 y2 

*1> *2’ I 3’ x 4 ’ ' 
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has a x distribution of one degree of freedom. 

Once a series of Eg’s, or resonance energies, 
from the appropriate x^ distribution; the value of 
then 


have been generated, we can sample 

XL 

r n associated with the n Eq is 


r n<v = E 0 / 2 r^0 < s -' rave > 


( 9 ) 


where Vq = 1. 


Program DAISY 

The program DAISY (ref. 8) is given the average level spacing for a resonance 
chain and the average reduced scattering width. It then generates many pseudo- 
resonance chains or sets of resonance parameters (as many as 9999 chains of at most 
5000 resonances) with level spacings chosen from Wigner’s level spacing distribution 
and scattering widths chosen from the distribution of Porter and Thomas as described 
in the preceding section. The code then rejects those chains whose averages are beyond 
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the limits associated with a preselected probability level. This allows us to select a 
confidence probability level, say 95 percent, and to reject a chain if it does not approx- 
imate the average level spacing and/or the average scattering width better than 95 per- 
cent of all chains generated. This process limits the allowable variation in the average 
cross sections computed from such selected chains. The code then chooses, from these 

selected chains, that particular chain whose individual parameters best fit the theoreti- 

o 

cal distribution functions. Either the Kolmogorov-Smirnoff or \ tests for goodness of 
fit are used for such a selection of this best chain. Thus, the chain finally chosen by 
the program DAISY shall have the correct average values, and the distribution of the 
chain's individual parameters will most closely resemble their theoretical distribution 
functions. 


CROSS SECTIONS AND RESONANCE PARAMETERS FOR TANTALUM 

In order to use the method of pseudo-resonance parameter chains, we must have 
values of the average reduced neutron width T the average level spacing D, the po- 
tential scattering cross section cr j., and the average capture width r^. Certain of 
these parameters will be varied parametrically as part of the analysis. However, ac- 
curate values of the averages and their associated experimental errors are necessary 
as they are used as constraints on the final results of the analysis. These parameters 
may be obtained from measurements in the resolved resonance region. The following 
section gives the source of the parameters used for tantalum, along with a brief evalua- 
tion of their accuracy. The average values and their associated accuracy chosen for the 
present analysis are shown. 


Potential Scattering Cross Section 

Two experimental techniques which allow the determination of the potential scatter- 
ing cross section are the bright-line technique and the average-cross-section method. 

In the "bright- line" technique, the total cross section is measured by a high-resolution 
crystal spectrometer near a prominent resonance; for tantalum the 4. 28-eV resonance 
was used. The measured values are fitted by a resolution- and Doppler- corrected theo- 
retical cross-section curve. The parameters of the best-fitting theoretical curve are 
then reported as the resonance parameters for that particular resonance. Since both a 
resonance peak and a resonance wing analysis are performed, the potential cross sec- 
tion may also be extracted. In 1956, Wood, using the (BNL) crystal spectrometer, re- 
ported high- resolution measurements on four sample thicknesses of tantalum (ref. 9). 


6 


TABLE I. - COMPARISON OF WOOD (1956) RESONANCE PARAMETERS OF 4.28-eV RESONANCE 

IN 181 Ta WITH LATER MEASUREMENTS 


Parameter 

Source 

Wood (1956), 
ref. 9 

Evans (1958), 
ref. 10 

Doii'nitsyn (1961), 
ref. 11 
(a) 

Stolovy (1963), 
ref. 12 

Spin quantum number, J 

4 

4 

(4) 

4 

Statistical factor, g 

9/16 

9/16 

(9/16) 

9/16 

Scattering width, T n , me V 

3. 9±0. 6 

3. 6±0 1 

(3. 6±0. 1) 

— 

Capture width, r^, meV 

49±6 

62±2. 6 

51±5 

— 

Total width, T , meV 

53 ±5 

66 ±2. 7 

56±5 

— 

Resonance energy, E, eV 

4. 282±0. 008 

4. 28±0. 01 

4. 28±0. 01 

— 

Potential cross section, barns 

7. 0±1 0 



— 


a Only 2gr n reported, thus the parameters in parentheses were derived by assuming J = 4. 


the parameters extracted by Wood are shown in table I. As can be seen, later experi- 
ments (refs. 10 to 12) have confirmed the spin assignment J of the 4. 28-eV level. It 
should also be noted that later experiments have obtained somewhat higher values of the 
total width r , indicating the experimental value reported of CTp^ = 7. 0±1. 0 barns may 
be high. 

A second technique for the measurement of the potential cross section is the so- 
called "average- cross- section method. " If the transmission of neutrons averaged over 
many resonances, in the unresolved region, is plotted as a function of sample thickness, 
the slope for thin samples is proportional to the s-wave strength function. However, 
for thick samples the average transmission is 


T 





( 10 ) 


if it is assumed that the statistical factor g equals 1/2, r n = r, and resonance- 
resonance interference terms cancel. In equation (10), T p , is the transmission corre- 
sponding to a potential cross section a pot only, D is the average spacing for levels of 
a single spin state, nt is the sample thickness in centimeters, and ± is the neutron 
wavelength divided by 2n. Taking the logarithm, differentiating with respect to nt and 
neglecting higher order terms yields equation (11). 


— (In T) 
d nt 


°pot + 



( 11 ) 
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Thus, for thick samples the slope of In T as a function of nt is asymptotic to 

In an analysis of tantalum neutron transmission experiments, Seth, Hughes, Zim- 
merman, and Garth (ref. 13) in 1958 selected four "typical" resonances and computed 
their cross section, including interference effects and Doppler broadening. Then aver- 
age transmissions were computed for various sample thicknesses, and reasonable 
agreement with experiment was obtained for a = 8 . 5±1. 0 barns. 

Using the same technique, Divadeenam (ref. 14) in 1967 reported an extremely 
small value of the potential radius, which would imply a potential cross section of 
4. 80 barns. However, a very high value of the s-wave strength function Sq was also 
reported. The Divadeenam value of Sq for tantalum reported was 3. Ixl 0~4 (ref. 14 ). 
This may be compared with the Sp-value reported by Desjardins, Rosen, Havens, and 
Rainwater of 1. 84x10” 4 obtained from the individual measurement of 62 separated reso- 
nances (ref. 15). Thus, the Divadeenam value for the potential cross section may be 
somewhat suspect. 

The value for the potential cross section chosen for the present analysis is 8 . 5 
(+1.0, -2.5) barns. The negative value of the error is chosen to reflect the lower value 
of Wood (ref. 9). 


Average Level Spacing 

The average level spacing D of tantalum is determined by measuring the energies 
of individual resonances, subtracting the lowest resonance energy from the highest 
resonance energy measured, and dividing by one less than the number of resonances, or 


n - 1 

i-V. 

where E n is the resonance energy of the n n resonance measured from the first reso- 
nance energy Ej. A value of D for tantalum of 4. 5 eV was reported by Garg, Rain- 
water, and Havens in 1966 (refs. 16 and 17). These measurements extended from about 
30 eV to over 1300 eV and comprised over 300 neutron resonance energies. Figure 1 
shows the cumulative probability distribution of the level spacings measured experimen- 
tally compared to the theoretical cumulative probability distribution of two uncorrelated 
Wigner level spacing distributions. The cumulative probability distribution function 
F(s) of two Wigner single- species level spacing distributions which have been superim- 
posed on one another is given by the following formula: 
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0 .5 1.0 ' 1.5 2.0 2.5 3.0 3.5 4.0 


Normalized level spacing, D/D 

Figure 1. - Cumulative probability distribution of normalized level spacing 
compared to Wigner level spacing distribution for two spin species - for 
first 300 level spacings reported in reference 17. Number of level spac- 
ings, 300; ratio of average level spacings for two spin species, 1. 286; 
average level spacing of first 300 level spacings, 4.526 eV; energy Eq 
of last resonance.used in average and display, 1362. 2 eV. 



where k is the ratio of the average spacings of the two individual Wigner single-specie 
distributions (k = Dg/Dj) and s is the spacing. 
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If a (2J + 1) level density obtains, 7/16 of the levels should have J = 3 and 9/16 of the 
levels J = 4, where J is the spin of the compound nucleus. Thus, the appropriate 
value of k, or D 2 /Dp is 9/7 or 1. 286. As can be seen in figure 1 there is an excess 
of measured level spacings in the region between D and about 2. 5 D, and a deficit of 
measured level spacings in the region 0 to D. This indicates that if the Wigner model 
applies, many closely spaced resonances have not been resolved in the experiment of 
references 16 and 17. The fit of the experimental data to the theoretical curve is nearly 
the same for k = 1. The fit is somewhat improved for very large (or very small) k, 
when the theoretical curve approximates the single-species Wigner cumulative distribu- 
tion. However, in the low-energy region, from about 0 eV to about 300 eV, where few 
if any levels are missed experimentally, the fit of experimental data to the expected 
theoretical curve is quite good (see fig. 2 and appendix A); perhaps only three levels are 
missing. Using this low-energy, portion D is estimated to be 4. 3±0. 2 eV. 

Average Reduced Scattering Width 

The strength function Sq is generally defined as the ratio of average reduced scat- 



Figure 2. - Cumulative probability distribution of normalized level spacing 
compared to Wigner level spacing distribution for two spin species - for 
first (ft level spacings reported in reference 17. Number of level spacings, 
69; ratio of average level spacings for two spin species, L 286; average 
level spacing of first 69 level spacings, 4.344 eV; energy Eq of last 
resonance used in average and display, 304.0 eV. 
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Porter-Thomas distribution 
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Normalized reduced scattering widths, fjy l~^ 

Figure 3. - Cumulative probability distribution of normalized reduced scattering widths compared to Porter-Thomas distribution of reduced scattering widths - for 75 widths reported 
in reference 17. 



— 0 

tering width r for a particular J-state over the average level spacing for that partic- 

n — 0 0 

ular J-state. We assume that T n (J = 3) = T n (J = 4); and figure 3, since there is good 

agreement between the theoretical distribution and the experimentally measured distri- 
bution, indicates this to be a reasonable assumption for tantalum. If we further assume 
D(J = 3) = D(J = 4) = 2D, we may easily obtain the value of by the following equation: 

FjJ = 2 DS q (14) 

Two rather different means have been used to measure the s-wave strength function. 
One technique has been to measure the average energy-dependent behavior of the total 
cross section, to apply self- shielding corrections, and then to choose the s-wave 
strength function which best fits the resulting total cross-section curve. The other tech- 
nique has been to measure the individual r 's of resolved, low-energy resonances and 

n 0 0 
then to determine the strength function directly from the inferred values of T n or2gT n . 

Table n gives a summary of measurements of the s-wave strength function (refs. 14, 

15, and 18 to 21), along with the value of T*? obtained by use of equation (14). 

— 0 n 

A value of T near the mean of the values shown in table n has been chosen, and 

n 0 

an error has been chosen to cover the range of the experiments. A value of T n = 2. 0 

(+0. 9, -1. 2) meV has been chosen. 


Average Capture Width 

— o 

The average capture width is expected to be distributed as x with a large 
number of degrees of freedom and should be well approximated by the normal distribu- 
tion. The capture width has been measured for 26 resonances of tantalum. The weight- 
ed mean of these individual T^'s has been reported as 60 meV (ref. 15). However, 
only relatively large values of are measured since must be larger than or at 

least of the order of T to be measured. Hence, the distribution of measured r 's 

n 7 y 


TABLE II - MEASURED s-WAVE STRENGTH FUNCTIONS FOR TANTALUM AND ASSOCIATED SOURCE 



Source 


Simpson (1958), 
ref 18 

Desjardins (1960), 
ref 15 
(a) 

Seth (1964), 
ref 19 

Divadeenam (1967), 
ref 14 

Frike (1971), 
ref 20 

Dilg and Vonach (1971), 
ref 21 

s-Wave strength function 
Average reduced scattering 
width, r* meV 

(2 00*0 24)xl0~ 4 
1 72±0 20 

(1 84±0 34)xl0" 4 
1 58±0 29 

(2 50±0 50)xl0‘ 4 
1 58±0 29 

3 IxlO' 4 
2 67 

(2 2±1 2)xl0 4 
1 89±1 03 

2 6±0 5 
2 24±0 43 


3 Sq determined from individual resonances 
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Figure 4. - Distribution of measured capture widths of tantalum, along with normal distributions centered on 
mean and mode of measured distribution (refs. 15 and 17). Modal value is taken as more representative of 
mean of capture widths of tantalum (see text). Measured values are shown as histograms. 


is not the same as the expected distribution of T^'s but is biased toward higher values 
of r , figure 4. Thus, a measure of central tendency other than the mean should be 
chosen in order to eliminate some of the bias. Here the mode, or most probable value, 

has been chosen as a less biased estimator of the true value of r . A value of T = 

y y 

52±7. 5 meV is thus obtained. Normal approximations to the distribution of ob- 
tained by using both the mean and mode as the measure of central tendency are shown 
in figure 4. 


Nominal Values of Parameters 

The nominal values chosen for the average resonance parameters for tantalum are 
as follows: 

(1) Potential cross section, Up^, 8. 5 (+1. 0, -2. 5) barns 

(2) Average level spacing, D, 4. 3±0. 2 eV 

(3) Average reduced scattering width, T 2.0 (+0. 9, -1. 2) meV 

(4) Average capture width, r , 52±7. 5 meV 

These average values are used as initial values for the analysis and their associated 
errors are used as constraints of the analysis. 
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DESCRIPTION OF EXPERIMENTS AND ANALYSIS 


Description of Experiments 

Brief descriptions of the two experiments on tantalum at 2-keV neutron energy fol- 
low. If the reader is interested in more experimental detail, the original reports should 
be consulted. 

Rensselaer Polytechnic Institute (RPI) experiment - 1971 . - A series of neutron 
transmission and self-indication experiments were performed on tantalum below 100 keV 
(transmission results in fig. 5) to study the effects of resonance self-shielding in the 
unresolved neutron energy range (ref. 5). Only those results about 2-keV neutron 
energy, shown in table HI, are considered in this report. 

Room -temperature (295 K), high-temperature (1073 K), and low-temperature (77 K) 
transmission and self-indication experiments were carried out with sample thicknesses 
of 0. 0278 and 0. 0549 atoms per barn. These samples provided significant temperature 
variations of the effective cross section and correspond to nominal tantalum thicknesses 
of 0. 5 and 1. 0 centimeters. 

Likewise, self-indication ratios (SIR) were measured as a function of neutron energy 
and the temperature of the transmission sample. The self-indication experiment con- 
sists of a transmission sample of tantalum, which may be heated or cooled, and a detec- 
tor sample of tantalum which is inside the RPI 1. 25-meter liquid scintillator capture 


Temperature, 



Figure 5 - Average transmission over energy region 300 eV 
to 100 keV for two samples of tantalum. (Data from ref. 5, 
RPI 1971 ) 
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TABLE IB. - RPI 1971 RESULTS AT APPROXIMATELY 2keV FOR 


TANTALUM IN TRANSMISSION AND SELF-INDICATION 
(1825.0 TO 2017.0 eV) 

[Errors shown are three-sigma limits.] 


Slab thickness, 
cm 

Effective average total cross section, c^, barns, at - 

77 K 

295 K 

1073 K 

Transmission 

0 5 

15. 25±1 . 14 

15. 67±0. 83 

16. 59±1.08 

1 . 0 

13. 13±0. 52 

13. 90±0. 41 

15. 09±0. 58 

Self -indication 

0.5 

28. 90±2.74 

29. 98±1. 23 

30. 72±2. 20 

1.0 

25. 31±1. 03 

26. 41±0. 65 

26. 06±1.27 


detector. The SIR'S are then the ratio of captures in the detector sample with the trans- 
mission sample in the neutron beam to those with the transmission sample out of the 
beam. For these experiments the detector sample, which is down beam from the trans- 
mission sample, was 0. 00147 atom per barn thick, which corresponds to a thickness of 
0. 0266 centimeter and was maintained at room temperature throughout. 

Thus, for each energy region and sample thickness there are six independent 
measurements: three temperatures in transmission, and three temperatures in self- 
indication. 

Idaho Nuclear Corporation (INC) experiment - 1968 . - A series of neutron trans- 
mission experiments designed to determine the variation in total cross section as a 
function of temperature at 2 keV were performed at the Idaho Nuclear Corporation's 
Materials Testing Reactor in 1968 (refs. 6 and 7). A scandium filter was used as a 
monochromator for 2-keV neutrons. The minimum value of scandium's total neutron 
cross section at 2 keV allows its use as an energy filter for neutrons. A 84- centimeter 
(33-in. ) thick scandium filter was placed in the INC Materials Testing Reactor's fast 
chopper beam, and the transmitted neutron spectrum was measured. Part of the 
scandium -filtered spectrum is shown in figure 6. 

The 2-keV beam should make it possible to observe the temperature dependence of 
neutron cross sections caused by the Doppler effect over the energy width of the beam. 
The resolution of the beam is broad enough to contain a sizeable number of resonances 
(~200 levels), thus providing a measure of average behavior. 

The neutron transmission of a 0. 508-centimeter (0. 200-in. ) thick tantalum sample 
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was measured as a function of temperature. The percentage increase in total cross 
section as a function of temperature as measured in the INC experiment is shown in fig- 
ure 7. 

COMPUTER SIMULATION OF RPI 1971 EXPERIMENTS 

The results of the RPI 1971 experiments have been reduced to effective average 
total cross sections in the following manner: An effective average total cross section 
may be directly calculated from the average transmission T by 

= (^j In (T) (15) 

where n is the number density of nuclei per cubic centimeter and t is the transmission 
sample thickness in centimeters. Another effective average total cross section may be 
defined from the average self-indication ratio S by the equation 

(i6> 

where t is the thickness of the shielding sample in centimeters. It can be shown that 
the effective total cross sections are related to the true average total cross section a 
by the inequality 

(c t ) > a > ff t (17) 

s 

see appendix B. 

Given a set of DAISY-generated spin- dependent resonance parameters over the 
energy range of average cross sections, we may compute continuous cross sections by 
using a Doppler-broadening computer code (ref. 22). Since tantalum has a substantial 
Debye temperature (refs. 23 and 24), 0 D = 231 K, we must use an effective temperature 
rather than the ’’true" thermodynamic temperature in the Doppler-broadening computa- 
tion of the cross sections. The method of Lamb (refs. 25 and 26) has been used to gen- 
erate the effective temperature. Tantalum is assumed to be a Debye solid. The follow- 
ing formula was solved by numerical quadrature for the effective temperature T* : 

(18) 
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TABLE IV. - EFFECTIVE TEMPERATURE T* (K) OF 
TANTALUM AS FUNCTION OF THERMODYNAMIC 
TEMPERATURE T(K) FOR DEBYE 
TEMPERATURE OF 231 K 


T 

T* 

T 

T* 

1 

86.6 

400 

406 6 

50 

93 5 

500 

505 3 

100 

125. 1 

600 

604.4 

150 

167.3 

700 

703.8 

200 

213.1 

800 

803.3 

250 

260.6 

900 

903.0 

300 

308,8 

1000 

1002.7 


where T is the thermodynamic temperature. Table IV displays effective temperatures 
of tantalum for selected thermodynamic temperatures. It can be seen that reducing the 
thermodynamic temperature of tantalum below that of liquid nitrogen would not reduce 
the effective temperature of tantalum significantly. The effective temperatures used in 
the RPI 1971 computer simulations were 110. 5, 300. 1, and 1075. 5 K, for the liquid- 
nitrogen-temperature, room-temperature, and high-temperature samples, respectively. 

Two separate chains are generated for both s-wave J-states of tantalum. This 
procedure is followed to preclude the effect of merged cross-section sets (ref. 27). 

For these two separate chains, r ®(J = 3) is made equal to r^(J = 4); and likewise with 
D(J = 3) and D(J = 4), which are set equal to 8. 6 eV. Once these continuous cross sec- 
tions have been calculated, we may use these cross sections in a computer simulation 
of the RPI 1971 experiments. 

The average transmission of a 1/E neutron beam through a slab of material in an 
energy region from Ej ow to E^^ may be computed by the following formula: 

N 

_ X(i)H nt0<Ei ® 

T (19) 

N 

Eft) 

i=l 

where E i = E 1qw + [UE^^ - E 1qw )/n], a(E i ) is the value of a t (computed) at Ej and 
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N is the number of cross sections tabulated in the interval [E low , E highJ- This ex- 
pression, formula (19), is essentially exact for the RPI 1971 experiments, as any 
neutron that is scattered by the transmission sample is removed from the beam and not 
counted. The transmission so computed is transformed to a by equation (15) and 
can be compared directly with the measured a^'s. 

The average self- indication ratio S of a 1/E neutron beam through a slab of thick- 
ness t onto a thin self-indication sample of thickness t* may be approximated by the 
formula 


N 


s 


/ (— )expf-na(E.)t]-^ 

/ y \eJ l 1 J cr 


V E i) 


(EJ 


|l - exp[na t (E.)t*]J 



( 20 ) 


where a^(E^)is the total cross section of the self-indication sample evaluated at E^, 
o^(Ej) is the capture cross section of the self- indication sample evaluated at E^ and 
all other symbols have the same meaning as in equation (19). Equation (20) does not in- 
clude corrections for multiple scattering in the self-indication detector sample and is 
thus an approximate relation. The approximation involved in equation (20) has been 
checked by a single scattering model and by Monte Carlo simulation of multiple scatter- 
ing in the detector sample. Equation (20) is found to be quite accurate, within 0. 5 per- 
cent in calculating the self-indication ratio. 

The value of S obtained from formula (20) may be converted into a (a. ) by equa- 

_ s s _ 

tion (16) and can be compared to the (cx^) ' s resulting from measurements of S. Recall 

s 

from inequality (17) that these two values, computed from transmission and(c^) 

computed from the self-indication ratio, should bracket the true average total cross 
section. 

Using formulas (19) and (20) and DAISY- generated resonance parameters, the RPI 
1971 experiments may be accurately simulated on a computer. 


Temperature-Dependent Analysis of RPI 1971 Experiment 

In this section the analysis of the RPI 1971 temperature-dependent experiments at 
approximately 2 keV is detailed. Figure 5 illustrates the transmission data obtained in 
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the RPI 1971 experiments. The precision of the experimental results is such that struc- 
ture caused by changes in the local strength function is readily discernible in the unre- 
solved region. It is clear that if we wish to model the detailed cross-section behavior 
actually measured, we must attempt to extract the local strength functions and an ana- 
logous set of resonance parameters. The region about 2 keV has been chosen for this 
analysis as it appears to be a region of particularly high transmission and consequently 
low strength, thus providing a severe test of the technique. The 2-keV region has also 
been chosen since there exists other temperature -dependent neutron transmission data 
for tantalum at this energy (refs. 6 and 7). 

Since the energy region chosen for pseudo-resonance analysis is rather narrow 
(192. 0 eV at 4. 3 eV per level, or 45 levels), we should not expect that the local average 
of r ® over these 45 levels will be exactly In fact, if we assume several values 

of T ® from the range allowed for the global average of r®, the probability table 
(table V) indicates the probability levels of various results. For example, if the global 
value of T"® is 1.5 meV, the probability of a local average of r® for 45 levels less 
than 1. 4 meV is 40 percent. As can be seen from the probability table, the range of 


TABLE V. - PROBABILITY TABLE OF LOCAL AVERAGE REDUCED 
SCATTERING WIDTH (local) FOR DIFFERENT VALUES OF 


AVERAGE REDUCED SCATTERING WIDTH T 0 

n 


Probability of local average of 
45 T^'s being greater than 

n Q 

or equal to T” (local) 

r£ (local), meV, if rj 

is - 

1. 0 meV 

1. 5 meV 

2. 0 meV 

2 5 meV 

0.001 

0.471 

0.706 

0.941 

1 177 

.005 

.539 

809 

1 079 

1 349 

.010 

.575 

863 

1 150 

1 438 

.050 

.680 

1 020 

1.360 

1 701 

100 

.741 

1 112 

1.482 

1 853 

200 

.820 

1 230 

1.640 

2 050 

.300 

880 

1 320 

1.760 

2 200 

400 

.933 

1 400 

1.867 

2.333 

500 

.985 

1.478 

1.970 

2.463 

.600 

1.039 

1.559 

2.078 

2.598 

700 

1.099 

1.648 

2 198 

2 747 

800 

1. 172 

1 757 

2 343 

2 929 

900 

1.278 

1.917 

2.555 

3 194 

.950 

1.370 

2 055 

2 740 

3. 425 

990 

1 555 

2.332 

3 110 

3.888 

.995 

1 627 

2 440 

3.253 

4 067 

999 

1.781 

2.672 

3. 563 

4 453 
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r ® (local) is rather large even for a region containing 45 levels. Thus, a method has 

n — 0 
been adopted in this study which will determine the local value of r . This is done by 

0 n 

seeking regions wherein T n (local) and c; ^ provide temperature-dependent effective 

total cross sections which are compatible with the RPI 1971 experimental results and 

minimize the difference between the local average and the measured averages. 

The computer simulation can compute quantities similar to those measured, that is, 
the average transmission and the average self -indication ratio. We can present the 
simulation results, as well as the measurements in terms of average effective total 
cross sections, by the use of equations (15) and (16). The average effective total cross 
sections generated by the computer simulation can be represented as a resonance- cross- 
section contribution plus the potential cross section. Several values of r jj (local) have 
been used in these computer simulations. The results of these simulations are shown 
in tables VI and VII. If we compare the nominal values shown in tables VI and VH with 
the measured values of table m, we realize that the use of global averages in this region 
predicts a Doppler effect of more than double that observed and a total cross section 
about 30 percent larger than that observed. Thus, it is clear that, to adequately model 
the cross-section behavior, local values must be used. 

Subtracting the potential contribution from the average effective total cross section 
leaves what will be called the resonance contribution to the average effective total cross 
section. This resonance contribution to the total cross section has been spline inter- 
polated and plotted as a function of the local value of T® for the three temperatures: 
cold, room, and high. In figures 8 and 9 the resonance portions of the average effective 


TABLE VI. - AVERAGE VALUES OF LOCAL PSEUDO-RESONANCE PARAMETERS AND RESULTANT EFFECTIVE TOTAL 
CROSS SECTIONS COMPUTED FOR 0 5-CENTIMETER-THICK TANTALUM SLAB 



Nominal 

Variations from nominal case 


2 

3 

4 

5 

6 

7 

Potential cross section, a pot , barns 

8.5 

8.5 

8 5 

8.5 

9 0 

9.0 

9.0 

9.5 

Average reduced scattering width, r eV 

002 

.0015 

.001 

. 00086 

. 00086 

. 00086 

0005 

.0007 

Average capture width, T^, eV 

.052 

.052 

.052 

.052 

052 

.045 

.052 

.052 

Effective average total cross section, a t , barns, at - 









77 K 

18.50 

16.75 

14.80 

14. 19 

14.69 

14.69 

12.88 

14.42 

295 K 

20 41 

18.31 

15.89 

15.12 

15 12 

15.61 

13.37 

15. 17 

1073 K 

23.37 

20.51 

17.25 

16.24 

16.75 

16.75 

13.88 

16.01 

Effective average total cross section defined from 









average self- indication ratio, (cj , barns, at - 









77 K 

44.57 

40.28 

34.77 

32.85 

33.34 

33.15 

27.05 

31.29 

295 K 

47.97 

42.64 

35.83 

33.51 

34.00 

33.70 

26.68 

31.51 

1073 K 

45.79 

39.63 

32.22 

29.83 

30 36 

30. 10 

23.23 

27.85 
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Figure 8. - Resonance portion of effective total cross section as determined from trans- 
mission through 0. 5-centimeter-thick slab of tantalum as function of local value of 
average reduced scattering width between 1825 0 and 2017.0 eV. Average level spacing, 
4.3 eV; average capture width, 0.052 eV. 



Figure 9. - Resonance portion of effective total cross section as determined from transmis- 
sion through 1 0-centimeter-thick slab of tantalum as function of local value of average 
reduced scattering width between 1825.0 and 2017.0 eV. Average level spacing, 4.3 eV; 
average capture width, 0.052 eV. 
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total cross section determined by transmission have been plotted as a function of 

rj? (local). In figures 10 and 11 the resonance portions of the average effective total 

n — 0 

cross section determined by self- indication have been plotted as a function of r n (local). 

It should be noted that the change of average effective total cross section with respect 

to a change in local strength is nonlinear. The smooth curves shown in the figures are 

spline approximations (ref. 28) obtained by using Akima's approximation (ref. 29) for 

the end-point slopes. 

In the same way, the measured, average effective total cross sections are also made 
up of a resonance portion and a potential cross section. Given these experimentally 
determined average effective total cross sections, we may subtract varying amounts of 
potential cross section and plot the remaining resonance portion and its error along the 
calculated resonance-cross-section curves of figures 8 to 11. For example, figure 12 
shows the average effective resonance portion of the total cross section as determined 
by transmission for the 0. 5- centimeter slab of tantalum at room temperature as a func- 
tion of rjj (local). The experimentally measured value of the average effective total 
cross section, for the 0. 5-centimeter room -temperature sample in transmission, is 
15. 67±0. 83 barns. For a potential cross section of 8. 5 barns the remaining 7. 17±0. 83 
barns may be plotted on the curve. The limits of the rectangular region of figure 12 
indicate the region of local r^'s which are compatible with a potential cross section 
of 8. 5 barns and a measured by transmission through a 0. 5- centimeter slab of 
tantalum at room temperature of 15. 67±0. 83 barns. The same process is repeated for 
all temperature measurements. If a region of consistent local T^'s is found with 



Figure 12. - Region of compatibility for 0. 5-centimeter-thick slab of tantalum at room temperature and potential 
cross section of 8. 5 barns. 
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^ Average local value of reduced scattering width, n* (local), eV 

Average local value of reduced scattering width, (local), eV 





respect to all temperature measurements of a sample and a potential cross section, 
this region may be plotted on a potential-cross-section and r® (local) chart, see fig- 
ures 13 to 16. These provide us with areas in which the pseudo-resonance parameters 
chosen by DAISY match the experimental results of the RPI 1971 experiment. 

It was found (cf. figs. 13 to 16) that the thick (1. 0 cm) sample of tantalum in trans- 
mission provided the most stringent region of compatible parameters. The thin (0. 5 cm) 
sample of tantalum in self-indication showed no region of compatibility with reasonable 
variations of the resonance parameters. However, it should be recognized that the self- 
indication regions of compatibility did not appear to be particularly sensitive to changes 
of the local average resonance parameters. The regions of compatibility for self- 



Figure 15. - Self-indication partially compatible regions for 0. 5-centimeter-thick tantalum slab Regions en- 
closed by curves will produce resonance cross sections which are compatible, as labeled, with RPI 1971 experi- 
ments on 0. 5-centimeter slab between 1825.0 and 2017.0 eV. 



Figure 16. - Self-indication compatible region for 1.0-centimeter-thick tantalum slab. Region enclosed by curves 
will produce resonance cross sections which are compatible with RPI 1971 experiments on 1.0-centimeter slab 
between 1825.0 and 2017.0 eV. 
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indication experiments are bands which are nearly horizontal. The regions of compati- 
bility for the transmission experiments are lunes of relatively small area. Thus, the 
transmission data have been given more weight in the choice of final compatible param- 
eters than the self-indication data. 

It may be seen from observation of figure 5 that the behavior of tantalum’s cross 
section is not smooth even with the 10 percent energy bins chosen. Also at 2 keV the 
local strength must be rather low as the transmission is enhanced. This is in agree- 
ment with the findings of this report and is caused by the local clustering or sparsity of 
resonances and/or the local strength of the resonances in these energy bins. Since 
there appear to be significant differences in the local strength functions, the present 
analysis demonstrates a technique for the extraction of the local strength function from 
temperature-dependent transmission and self -indication measurements. The technique 
demonstrated may be used to provide sets of pseudo- resonances whose individual param- 
eters obey the appropriate statistical distributions and whose average behavior simulates 
the transmission and self- indication measurements within experimental error. 

Tables VI and VII give the average values of the local pseudo-resonance parameters 
and the results of the computer simulations which use cross sections generated from 
these parameters. These tables allow us to ascertain the sensitivity of the effective 
total cross sections to changes in the pseudo-resonance parameters. It can be clearly 
seen from a comparison of variations 3 and 4 of the tables that, at least in this range of 
average parametric values, the changes in the effective total cross sections are directly 
proportional to the changes in the potential cross section. That is, a change of 
+0. 50 barn in the potential cross section produces almost a uniform increase of 
0. 50 barn in the effective total cross sections inferred from transmission and self- 
indication simulations. 

It may be observed by the comparison of variations 4 and 5 that the change in the 
average value of from 0. 052 eV to 0. 045 eV leaves the average effective total cross 
section inferred from transmission virtually unchanged. This change does alter the 
value of the average effective total cross section inferred from self-indication. This is 
caused by the changes produced in the capture cross section. The changes produced by 
altering T r are of second order compared to the changes produced by the local strength 
and the potential cross section. We can thus "fine tune" the final cross-section set by 
altering the average value of the capture width. 

As the self- indication experiments were subject to larger errors than the transmis- 
sion experiments, the values of a pQt (consistent) and (local consistent) were chosen 
to fit all transmission experiments and to fit most of the self-indication experiments. 

The values chosen to fit the RPI 1971 experiment in the energy region from 1825. 0 to 
2017. 0 eV were cr pot = 9.5 barns and T"® (local) = 7. 0x10"^ eV. These values, along 
with a r" of 52 meV, fitted within the error bars of all transmission experiments and 
four out of six self- indication experiments, figures 17 and 18. The result of a recent 
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Temperature of transmission sample, K 

Figure 17 - Comparison of measured and calculated average effective 
total cross sections, from transmission and self-indication ratio 
experiments, as function of temperature for 0.5'centimeter-thick 
tantalum slab. 


high- re solution transmission experiment is also shown in figure 17 (private communica- 
tion of 1968 data from F. J. Rahn and H. Camarda of Columbia University). 

A list of the DAISY-generated compatible resonance parameters for the energy re- 
gion 1825. 0 to 2017. 0 eV is shown in table VUI. This table is then the final result of the 
analysis of the RPI 1971 experiment. These parameters, while not unique, will ade- 
quately model the temperature -dependent behavior of tantalum at 2 keV and are suitable 
for use in unresolved resonance reactor calculations. 


TABLE VIII - RESONANCE PARAMETERS USED FOR COMPUTED 

CROSS SECTIONS SHOWN IN FIGURES 17 AND 18 

[Seventh variation used in RPI 1971 experiment (ref. 5), average 
reduced scattering width, r®, 0 70 meV, capture width, r^, 

52 meV; potential cross section, a 9. 5 barns] 


(a) Average level spacing, D(J = 4), 8 6 eV, 
statistical factor, g, 0.5625 


Resonance energy, 
E, 
eV 

Scattering width, 

r n’ 

eV 

1825 0000 

0.02990399 

1825 3337 

01448790 

1833 4493 

00087920 

1841. 6215 

. 04647899 

1851.3175 

00289850 

1866.7056 

.09095930 

1877.8835 

. 17080720 

1890.8943 

. 01664480 

1900.0273 

.07454510 

1905.3980 

.01845030 

1912 4215 

. 00064009 

1916 0287 

. 00000240 

1928.3157 

00040600 

1944.2893 

.00453179 

1949 8231 

.00001340 

1954. 1996 

00867499 

1962 6062 

00687890 

1967 5294 

04249540 

1977.7147 

.00701980 

1992.3089 

.07792800 

2009.2010 

.05289979 

2015 6339 

.02148139 

2019 5734 

03189699 

2022 6075 

. 01492579 


(b) Average level spacing, D(J = 3), 8 6 eV, 
statistical factor, g, 0.4375 


Resonance energy, 
E, 
eV 

Scattering width, 

r n- ! 

eV 

1825.0000 

0 02990399 

1832.6116 

.01888040 

1838 6334 

05318150 

1850. 6213 

.00016689 

1860.2224 

01661870 

1875.3686 

.02331850 

1892 0088 

02956620 

1900. 6927 

10772499 

1903 4424 

. 00773039 

1916 9520 

.05056860 

1922.2507 

.06293280 

1942.2597 

. 00004520 

1946.6656 

. 15881540 

1955.8728 

.01033830 

1959 7644 

06480910 

1964. 1843 

.00071040 ! 

1965. 1088 

.04471409 

1971 5433 

.03200819 | 

1982. 1331 

. 00262450 

1989.3610 

.00197519 

1997.3739 

.00122420 

2011.6603 

. 00064850 

2016.4285 

.00745690 

2022 6306 

.00777750 
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Temperature-Dependent Analysis of INC 1968 Experiment 


In reference 7 an attempt was made to fit the temperature-dependent change in the 
total effective cross section of tantalum at approximately 2 keV by a pseudo-resonance 
technique. The change in the total cross section as a function of temperature for several 
pseudo-resonance chains and for the experimental data is shown in figure 19. It was 
concluded in reference 7 that the pseudo-resonance curves computed, labeled A, B, C, 

D, and F, did not adequately describe the experimental data. Why the pseudo-resonance 
approach failed in this test poses a serious problem. 

Recall that the INC sample was 0. 508 centimeter (0. 200 in. ) thick; this allows us 
to compare the INC experiment directly with the RPI 0. 500-centimeter-thick sample. 
Since the scandium -filtered beam is considerably wider than the RPI energy resolution, 
the scandium beam shown in figure 6 has been broken into nine energy groups and folded 
over the RPI transmission data of figure 5. Two energy- averaged total effective cross 
sections were obtained over the scandium-filtered beam, one for room-temperature 
~o± = 18. 272 barns and one for the high-temperature (800° C) a. = 19. 662 barns. Thus, 
the measured scandium -filtered beam folded over the RPI-measured effective cross sec- 
tions gave a ratio of hot to room -temperature total cross sections of 1. 076 at 800° C. 



Figure 19. - Percentage change in average effective total cross section 
as function of temperature. Scandium -filtered beam of 2-keV neu- 
trons; average level spacing, 4.3; ratio of average reduced scattering 
width to average level spacing, 1 8x10 . Curves A, B, C, 0, and F 
represent five different sets of mock data that have been Ooppje r 
broadened as a function of temperature Point labeled RPI *%a 
represents the RPI effective cross sections folded over the scandium- 
filtered spectrum. (Data from ref. 7.) 
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This point has been plotted on figure 19. It is seen that this point is in good agreement 
with the "A-set" of pseudo cross sections generated in reference 7. Such agreement 
suggests that perhaps the experimental values of INC 1968 are subject to some source of 
experimental bias toward larger changes in the average effective cross section with 
respect to changes in temperature. This was, in fact, the case as the collimator into 
which the scandium was inserted had an intolerable amount of leakage through duct 
streaming (ref. 30). Thus, the tantalum sample was not only exposed to a 2-keV beam 
but also to a reactor- like beam of duct streaming neutrons which would serve to accen- 
tuate the Doppler effect measured. This deficiency has been noted and the scandium- 
filter - collimator system has been redesigned (ref. 31). However, the experiment has 
not been conducted with the new 107- centimeter (42-in. ) scandium-filter - collimator 
system. 


CONCLUSIONS 

A temperature-dependent pseudo-resonance parameter analysis has been performed 
on tantalum at 2-keV neutron energy. An earlier temperature -dependent analysis in 
this energy region has been discussed. The following conclusions were drawn: 

1. In portions of the unresolved neutron resonance region, where experimental resolu- 
tion is adequate to encompass few resonances «100) in a measurement the data permits 
the evaluation of the local strength function. Therefore, the local behavior of the meas- 
ured neutron cross section may be modeled by the use of appropriately chosen pseudo- 
resonance parameters. These parameters may then be used in reactor calculations with 

a measure of credibility similar to that for resolved resonance parameters. The cus- 
tomary approach of using pseudo-resonance parameters based on an average strength 
function rather than a local strength function can be erroneous. 

2. In the energy region 1825. 0 to 2017. 0 eV the temperature -dependent behavior of 
tantalum's effective total neutron cross section as determined by transmission has been 
adequately modeled by a specific set of pseudo-resonance parameters, using the tech- 
nique described. 

3. The technique illustrated in this report should be generally applicable to meas- 
urements in the unresolved neutron resonance region, provided that the inelastic and 
p-wave cross sections are either negligible or known. This technique provides a level 
of confidence that locally generated resonance parameters can model observed 
unresolved-cross-section behavior better than the naive sampling from global distribu- 
tions. 
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4. The results of this study suggest that a Doppler coefficient calculated by sampling 
from grand average statistical distributions over the entire unresolved region can be in 
error since local variations in the strength function are not taken into consideration. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 9, 1973, 

503-25. 
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APPENDIX A 


COMPUTER PROGRAM MOBEE 

Over 300 resonance energies have been measured for tantalum. To compare the 
theoretical level spacing distribution with the distribution of the experimentally deter- 
mined values, we could plot the cumulative distribution of all 305 level spacings and 
compare it to the theoretical level spacing for two uncorrelated spin species. However, 
we would find that the agreement is better with the single- species Wigner level spacing 
distribution than with the two- species distribution (figs. 1 and 20). This indicates that if 
the Wigner model applies, some closely spaced resonances have not been resolved in the 
region between 0 and 1397. 9 eV. We could laboriously plot the cumulative distribution 
for selected energy intervals and attempt to find the region in which serious deviation 
between the theoretical distribution and the experimental distribution arises. A some- 
what different approach has been adopted in the MOBEE computer code. It produces a 
plot of the experimental cumulative level spacing distribution and the theoretical cumula- 
tive level spacing distribution for each new experimentally measured level spacing in a 
35-millimeter motion-picture format. 



Normalized level spacing, D / D 

Figure 20. - Cumulative probability distribution of normalized level spacing 
compared to Wigner level spacing distribution for a single spin species - 
for first 300 level spacings reported in reference 17. Number of level 
spacings, 300; average level spacing of first 300 level spacings, 4. 526 eV; 
energy Eq of last resonance used in average and display, 1362. 2 eV. 
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W1GNER SINGLE SPECIES 



representative frames (in order) of output from computer code MOBEE. 






MOBEE, by the use of computerized microfilm plotting (ref. 32), produces a 35- 
millimeter motion picture of the experimental distribution's behavior relative to the 
theoretical distribution. By the use of a film editor, the viewer can easily see the point 
at which the experimental distribution begins to deviate significantly from the theoretical 
distribution. Sufficient fiducial data is listed on the film frame such that the D can be 
easily determined at that point. 

A series of resonance cards, containing at least the value of Eq for a resonance, 
are read in ascending order by MOBEE. At this point, several options are available to 
the code's user. He may specify an average level spacing D and let all distributions 
be plotted with this value as the mean value, or he may allow the code to compute aver- 
ages. The experimental cumulative distribution of level spacings can be displayed with 
no theoretical distribution, with a theoretical single -species cumulative distribution, or 
with a theoretical double-species cumulative distribution shown. If the user chooses the 
double-spin-species option, he must supply the ratio of D£/Dp of the two spin states. 

For illustrative purposes, a few frames from the MOBEE output are shown in fig- 
ure 21. 

The MOBEE code has been extremely useful in allowing the determination of the 
energy region in which closely spaced levels begin to be missed or not resolved. 
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APPENDIX B 


COMPARISON OF AVERAGE EFFECTIVE TOTAL CROSS SECTIONS FROM 
TRANSMISSION AND SELF-INDICATION WITH 
AVERAGE TOTAL CROSS SECTION 

The RPI 1971 experiments measured the average effective total cross section both 
from transmission and from self- indication. This appendix will indicate the relation 
of these quantities to the "true" average total cross section. This derivation follows 
that of Block and Byoun (ref. 32). 

Since the transmission sample is located in a neutron beam, the average transmis- 
sion (T) in the energy interval Ej to Eg is 


<T> 



(Bl) 


where n^ is the number of nuclei per cubic centimeter in the transmission sample, t 
is the thickness of the transmission sample in centimeters, and cr(E) is the total cross 
section as a function of energy. This may be rewritten as 


<T> 



^-n 1 t[o(E)+<o)-(o)] 


dE 


E 


2 


- E 


1 


(B2) 


where (o) is the "true” average total cross section between energies E^ and Eg- By 
rearranging expression (B2) we have 


<T> 


= e 


/ 2 e" n 

n^t(o) 


- ni t[a(E)-<o>] 


dE 


E 2 ‘ E i 


(B3) 
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or if we use brackets to again denote an energy average between energies Ej and E2, 

(B 4 ) 


- ni t<o> / - ni t[a(E)-<o>]> 
( T) =e (e 


We may then define 


V = (e 


-n 1 t[a(E)-< a)]' 


(B 5 ) 


as a cross-section fluctuation correction. If 


a(E) - <o> < 


(n t t) 


(B6) 


then we may approximate 77. 


„ = / e - n ltC<’(EH°>]\ ; 1+ I n 2 t 2 var[cr(E)] 


(B 7 ) 


since var[a(E)] = < ct^(E)} - ( 0) As the effective average total cross section is calcu- 
lated directly from the average transmission, we have 


( °t> 

\ n iV 


(B8) 


for the average effective total cross section from transmission, or 

<a.> =<o> - -Lin 77 (B 9 ) 

n l t 

From approximation (B 7 ), we find 77 £. 1 ; thus, 

<a t > < (o> (BIO) 

Here then we have shown the relation of the average effective total cross section as 
determined by transmission ( a^) and the average total cross section ( 0) . 
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The average self-indication ratio S may be defined as follows: 


S = 


< 


- ni to(E)\ 
Y y (E)e 1 / 

<Y y (E)> 


(Bll) 


where Y y (E) is the capture yield at energy E for the self- indication capture sample. 
The capture sample has ng nuclei per cubic centimeter and is tg centimeters thick. 
Neglecting multiple scattering effects within the thin capture sample (ngtg is only 
0. 00147/barn) the capture yield may be approximated by 


. Or (E) 
Y. .(E) «=-*— 


a(E) L 


-n2t2cr(E ) 


(B12) 


where <x y (E) is the capture cross section at energy E. Since n2t2<7(E) is small in the 
unresolved region, we can expand the exponential and approximate 


-n 9 t 9 a(E) 

1 - e A = ngtg^E) 


(B13) 


Hence, 


tfjE) 

Y '(E) =-£ — n 2 t 2 a(E) 
y a(E) 1 * 

= n 2 t 2 a y (E) (B14) 

Combining equations (B12) and (B14) yields the following approximation for the self- 
indication ratio: 
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/ -n 1 ta(E)\ 

^ 2 t 2 a y (E)e J 


<n 2 t 2 ^ y (E)> 


K(E)e 


-njt^E)^ 


( °y(E)> 


[Oy(E)e 


-n 1 t[cr(E)+< o) -< o)] 


< O'y (E )) 


= e 


-"it <o> \y 


c.,(E)e 


- ni t[a(E)-(o)] 


<°y(E)> 


(B15) 


Here again we can define the right-most term of equation (B15) to be a fluctuation cor- 
rection 4; that is, 


i =■ 


y°y(E )e 


-n 1 t[(7(E)-<o)]' 


<°y( E) > 


S = e 4 


(B16) 


(B17) 


Assuming a small fluctuation in the total cross section compared to l/(njt), we can ex- 
pand the exponential as in (B13) 


<a (E)[a(E)-(o)]> 

4=1- n,t — r - 

<°y(E)> 


(B18) 


The term 4 can be related to rj, the fluctuation correction for transmission. 
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(B19) 



Now let us consider the expansion of £ shown in equation (B18); the value of 
(a^(E)[a(E) - <o)]> is positive because the cross-section fluctuations are caused by the 
resonances. The positive contribution to this term where the capture cross section is 
near the peak of a resonance is greater than the corresponding negative contribution 
when the capture cross section is small; that is, far from the resonance peak. Thus, 

£ should be smaller than 1. 
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= — In (smaller than 1) 

(njt) 



( o) + Positive term 



Finally, 

(( < ’t) s ) a <°> 2<c t > 


(B23) 


or the true average total cross section is bracketed by the effective average total cross 
section determined by self-indication on the high side and that determined by transmis- 
sion on the low side . 
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